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The phenomenon of charged-particle oscillation in DC voltage biased plane-parallel conductors is discussed. 
Traditionally accepted mechanism for explaining the oscillatory behavior of charged particles in such system 
attributes the phenomenon to a process of charge exchange, which takes place when charged-particle is in close 
proximity to one of the electrodes. A novel finding here reveals that for microscopic or smaller particles under 
special circumstances, charged-particle oscillation does not involve charge exchanges. Such system radiates and 
the frequency of emitted radiation is controlled by a DC voltage biased across the two plane-parallel electrodes. 
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cumstances, core-shell structured charged-particle oscillates 
under constant electric field. Such oscillatory behavior can- 
not be explained by the aforementioned traditional picture be- 
cause it does not involve any charge exchanges. The criterion 
for such oscillatory modes is given by 



I. INTRODUCTION 

The phenomenon of charged-particle oscillation in DC volt- 
age biased plane-parallel conductors (or electrodes) is well 
known. 1 " 3 Such phenomenon has been extensively studied 
over years, both theoretically and experimentally, due to 
its usefulness in variety of applications such as electrostatic 
thruster and nanoprinting, for instance, which require highly 
energetic charged-nanoparticles with very high speedir— 

Traditionally, the phenomenon of charged-particle oscil- 
lation subjected to a constant electric field has been at- 
tributed to a process of charge exchange, which takes place 
when charged-particle is in close proximity to one of the 
electrodes. 9 When charged-particle is placed between a DC 
voltage biased plane-parallel electrodes, it migrates to the 
electrode of opposite polarity. For electrically conducting par- 
ticles, charge exchange takes place near the point of contact 
with the electrode. This reverses the polarity of the charged- 
particle and the particle gets repulsed towards the opposite 
electrode. There, again, the charge exchange occurs and this 
process gets repeated, resulting in charged-particle oscillation 
between the electrodes. Such process is schematically illus- 
trated in Fig. [TJ 
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Figure 1: (Color online) Schematic of charged-particle oscillation 
due to repeated charge exchange. 

A novel finding in this work reveals that under special cir- 
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Here, £o is the permittivity of free space, g = 9.8 m • s is 
the gravitational constant, h is the gap between the two plane- 
parallel electrodes, Qj is the positive effective charge carried 
by the core-shell structured particle, m is the mass of the parti- 
cle, and K 3 is the dielectric constant for the space between the 
two electrodes. The core-shell structured particle, which is il- 
lustrated in Fig. |2] has a outer radius of b, the inner conductor 
core radius of a, and the insulator shell layer has a dielectric 
constant of K%. The term Zdm is the maximum value assumed 
by the parameter Zd in Fig. |2]as the particle executes an oscil- 
latory motion. The particle oscillates back and forth between 
Zd = and Zd = Zdm\ hence, physically, Zd,m represents the 
turning point where the core-shell structured charged-particle 
begins to move back towards the conductor plate located at 
Zd = 0. 

For a negatively charged core-shell structured particle, the 
oscillatory criterion is given by 
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The two criteria, Eqs. ((TJ and @, are not form wise identi- 
cal due to the fact that negatively charged particle oscillates 
near the lower conductor plate whereas the positively charged 
particle oscillates near the upper conductor plate. Here, the 
upper conductor plate, which is located at z c i = 0, is assumed 
to be held at a higher voltage than the lower conductor plate, 
which is located at Zd = h. The finding in this work reveals 
that for the configuration illustrated in Fig. [2] the positively 
charged core-shell structured particle only oscillates in the re- 
gion where < z c ( < h/2; and, the negatively charged core- 
shell structured particle oscillates only in the region where 
h/2 < Zd < h. 

This remarkable result, which I find it quite obvious but not 
trivial, is a direct consequence of solving electrostatic bound- 
ary value problem involving core-shell structured charged- 
particle subjected to DC voltage biased plane-parallel elec- 
trodes depicted in Fig. [2] The problem is analyzed by first 
solving the electric potential V3 in the region M3 . The solu- 
tion to V3 is obtained by solving Laplace equation with appro- 
priate electrostatic boundary conditions. Thereafter, induced 
charges on the surface of each plane-parallel electrodes are 
computed. Similarly, the effective charge carried by the core- 
shell structured particle is also computed. The dynamics of 
the charged-particle is formulated by considering Coulomb 
forces between the charged-particle and the induced charges 
on the surface of each plane-parallel electrodes. To generalize 
the problem to include all particle speed ranges, the dynamics 
is formulated in relativistic formalism. 

This article is organized as follow: (I) introduction, (II) out- 
line of results, (III) theory, (IV) concluding remarks, and (V) 
acknowledgments. Normally, the section of theory immedi- 
ately follows the introduction. However, as this work involves 
lengthy derivations, albeit straightforward, the reader is prone 
to lose the essence of what this paper tries to portray. For this 
reason, the section titled "outline of results" is placed immedi- 
ately after the introduction. In the outline of results, no details 
of derivations are provided. Instead, the essence of this article 
is briefly summarized there using only the results, which are 
rigorously derived in the theory section. 



II. OUTLINE OF RESULTS 

The essence of this article is to investigate the dynamics of 
charged-particle illustrated in Fig. |2a), where a conducting 
sphere coated with a shell of dielectric is placed in an other- 
wise constant electric field. The upper conductor plate elec- 
trode is held at a DC voltage of Vr and the lower conductor 
plate electrode is held at DC voltage of Vl, where Vi <Vr- 
The two plane-parallel conductor plate electrodes are sepa- 
rated by a gap of h. The conducting sphere coated with a shell 
of dielectric has a free charge densities of <7i and 02- where 
<7i is the free charge density on the surface of the conducting 
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Figure 2: (Color online) (a) Cross-sectional view of an ionized core- 
shell structured particle confined by the DC voltage biased plane- 
parallel conductors, (b) Spherical polar coordinate system showing 
spherical polar triplet (r,6,tj>) of a vector R in Euclidean three di- 
mensional space, R 3 . 



sphere and 02 is the free charge density on the outermost sur- 
face of the dielectric shell. The free charge on the surface of 
dielectric shell has been introduced purely for generalization 
of the problem. Without loss of generality, 02 can be set to 
zero in the final form of solution. 

That said, adopting the particle coordinate system illus- 
trated in Fig. |2fb), the electric potential in region M3 of Fig. 
12 a) is given by 
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where C is a constant, E p is the magnitude of DC electric field 
in the gap between two parallel plates in the absence of core- 
shell structured charged-particle, 

E P = U\V T -V L \); 
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and, y and v are defined as 
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Here, Eq is the electric permittivity of free space and K 2 and 
K3 are dielectric constants respectively for regions M 2 and M3 
illustrated in Fig. |2] With V3 (r, 9) and application of appro- 
priate electric boundary conditions to each conductor plates 
yields 
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where a, Mp = <T, Mf , (p , s) is the induced surface charge density 
on the surface of upper conductor plate, Ou p = Ou p (p,s) is the 
induced surface charge density on the surface of lower con- 
ductor plate, and p = \J x 2 +y 2 . 

The net force exerted on the core-shell structured charged- 
particle by induced charges on each surfaces of the conductor 
plates is given by F = Fi +F2, 
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where i = (1,2) , q\ = <Ji up , g 2 = <Ju p , £3 is the electric permit- 
tivity of the region M3, and R, is given by 

Ri = e A pi cos 0i + e v pi sin <j>x + e z s, 

R2 = e x p2 cos 02 + QyPi sm <h + e z (s — h) . 

In the limit the charged-particle becomes very small com- 
pared to the dimensions of parallel plates, which is the case 
for micro- or nano-sized charged-particle confined in a micro- 
scopically large, but macroscopically small parallel plates, the 
F, of Eq. (|5]l for i = (1,2) can be shown to become 
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where 

Qt = %7ia(b-a)(Ji VAn(a 2 0\ +b 2 a 2 ) . 

K 2 

When the gravitational effect is included, the force experi- 
enced by the core-shell structured charged-particle is 

F7 = Fi + F2 — t z mg 



or 
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where m is the mass of the particle, g — 9.8 m • s is the grav- 
ity constant, and the gravitational force has been assumed to 
be in the — e z direction. It can be shown that Qt is related to 
V by 

Qt =47re K-3V; 
and the force Ft may be re-expressed, for convenience, as 
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It is noticed that F7-, which is the net force exerted on the 
core-shell structured charged-particle illustrated in Fig. |2fa), 
is a one dimensional force that only depends on the relative 
length, s, measured between the particle and the surface of 
the upper parallel plate electrode. 

The dynamics of oscillating charged-particle is given by 



where c — 3 x 10 8 m • s _1 is the speed of light in vacuum. Us- 
ing the explicit expression for Fy, Eq. (0, it can be shown 
that 
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where e z has been dropped for convenience and the notations s 
and s respectively denote the first and second time derivatives, 
i.e., s = ds/dt and s = d 2 s/dt 2 . 

In terms of the Zd parameter illustrated in Fig. 0a), 

s = z d + b, s = z d , s = z d , 

which is the separation length between the upper electrode 
plate and the uppermost surface of the core-shell structured 
charged-particle, the previous nonlinear ordinary differential 
equation becomes 
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To solve and plot Eq. ( [Tol l, the core-shell structured particle 
in Fig. |2] has been chosen to be the aluminum nanoparticle, 
where the core is aluminum and the shell is aluminum oxide. 
The following parameter values have been assigned: 

K 2 = 6, KT 3 = 1, 

a = 1 .5 |im, h = 1 mm, 

b — a = 4nm, 
Vr = 8kV, V L = 0V, 
CTi =0.014C-m- 2 , (LL> 

(7 2 =0C-m- 2 , 
p m ,i = 2700kg -m- 3 , 
p,„, 2 = 3800kg -m- 3 , 

where p m \ and p m2 are mass densities of the aluminum core 
and the aluminum oxide, respectively. The thickness of alu- 
minum oxide layer has been set at 4nm, which is typical 
of aluminum nanoparticles. 10 Because aluminum oxide is an 
high-k dielectric material, i.e., k 2 ~ 6, the a 2 has been set to 
zero.— For an insulator, the value of a 2 is negligible compared 
to <7i . The mass of the core-shell structured particle has been 
computed as 



m = -na p mt i + -n(b -a ) p m , 2 , 
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where m c and m s represent the masses of the core and the 
shell, respectively. With these values assigned for each of the 
parameters, Eq. (fTOt is solved via Runge-Kutta method sub- 
jected to the following initial conditions, 



z d (0) = 0.25/; and z d (0) = 0, 



(13) 



which conditions are schematically illustrated in Fig. [3] 

In this work, it is assumed that core-shell structured 
charged-particle is already ionized. And, for the initial con- 
ditions, this ionized particle just happens to be at Zd — 0.25/z 
with no initial speed. Then, the question to be asked in this 
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Figure 3: (Color online) Initial position of the particle. There are 
no free surface charges on the insulating shell; and, hence, a 2 = 
OC-nT 2 . 



work is this; what happens to the dynamics of this charged 
particle afterwards? Before I present the result of dynamics, 
I shall discuss how such particle might be ionized experimen- 
tally. 

Szirmai demonstrated that a core-shell structured particle, 
i.e., an aluminum core surrounded by a thin layer of aluminum 
oxide, can be charged by exposing it to a static electric field of 
sufficiently large strength. At the onset of the field emission, 
electrons are stripped from the particle's conductive core and 
tunnel through a thin dielectric layer, leaving the particle net 
positively charged as a whole. 12 Because electrons are more 
strongly bonded in a dielectric, it takes significantly larger 
electric field to ionize insulators. Therefore, it is reasonable to 
assume no free surface charges on the insulating shell in Fig. 
|3] implying o 2 = 0C • m~ 2 there. 

The field emission process strongly depends on the geom- 
etry as well as on the orientation of host material^"— For 
instance, a spherical surface has a lower field emission thresh- 
old point than a flat surface; which is also the reason why a 
conductive needle emits electron better than a thick rod con- 
ductor. Physically, a sharp tip can be described by a surface 
with large curvature whereas the dull one is described by a 
surface with smaller curvature. By definition, the curvature of 
a circle of radius r is large for small r and is small for large 
r. Accordingly, spherical nanoparticles have very large curva- 
ture whereas macroscopic spherical particles have very small 
curvature. What is referred to as a flat surface is just a spe- 
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cial case in which the radius r of a sphere becomes infinite 
in extent. Since the conductive needle emits electrons better 
than a thick rod conductor, a needle has a lower field emis- 
sion threshold than a rod. This implies that smaller spheri- 
cal conductors, such as nanoparticles, have lower field emis- 
sion threshold than larger spherical conductors or flat surfaces. 
That explained, the field emission thresholds are schemati- 
cally summarized in Fig. |4]for a spherical conductor, spherical 
dielectric, and a conductor plate, wherein the following field 
emission thresholds are assumed: (a) E cs for the conducting 
sphere, (b) E^s for the dielectric sphere, and (c) E cp for the 
conductor plate in vacuum. Since the geometry of conductor 
plates is just an extension of infinitely large sphere, the finite 
sized conducting spheres have lower field emission threshold 
than large conductor plates; and, thus, E cp > E cs . Because 
electrons are very strongly bonded in a dielectric, exception- 
ally large electric field must be applied to strip an electron 
from a dielectric material. In general, it is much easier to 
strip an electron from a conductor plate than from a dielectric. 
Hence, it is reasonable to assume E c / S > E cp . In summary, the 
three field emission thresholds satisfy the inequality given by 

Eds > E cp > E cs > 0, 

where Ej s , E cp , and E cs are field emission thresholds for di- 
electric sphere, conductor plate, and conductor sphere, respec- 
tively. 

What if the conducting spherical particle is coated with a 
thin layer of a dielectric shell, just like the one illustrated in 
Fig. [3j would it be still possible to ionize such particle by a 
process of field emission? The answer to this is yes, of course. 
For instance, Konopsky et al. have experimentally measured 
field emissions from sharp silicon tips covered with thin di- 
electric calcium fluoride layers. 18 Iwasaki and Sudoh investi- 
gated electron tunneling through an aluminum oxide thin film 
on a nickel-aluminum metal composite. 19 And, Kurnosikov et 
al. have also investigated electron transport through alumina 
oxide tunnel junctions. 20 Nevertheless, the insulating dielec- 
tric shell on the surface of spherical conductor increases the 
minimum electric field required to ionize a particle. Then, to 
ionize a core-shell structured particle, it is crucial that two 
plane-parallel conductor plates are able to produce electric 
field large enough to ionize a particle without electric dis- 
charge setting in. When there is a field emission originating 
from one of the conductor plates, i.e., the lower conductor 
plate in Fig. [3] the positively ionized core-shell structured 
particle quickly neutralizes. To prevent this, the two plane- 
parallel conductor plates must be able to produce sufficiently 
large electric field to ionize the particle, but at the same time 
this electric field must not be large enough to initiate field 
emission from conductor plates themselves. 

Zouache and Lefort investigated the phenomenon of elec- 
tric discharge in plane-parallel conductor plates, such as the 
one illustrated in Fig. [3] but without particle inside^ In 
their configuration, two plates are separated by an empty 
space gap of one micron in length. They have tested vari- 
ous materials for the conductor plates. Among various ma- 
terials tried for conductor plates, one was prepared from a 
mixture of 60% silver and 40% nickel in its material com- 
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Figure 4: (Color online) Schematic of electron field emission thresh- 
old. When the magnitude E p = 1 1 Ep 1 1 of applied electric field is less 
than the field emission threshold, spherical particle is polarized in- 
side, but no electrons are field emitted and the particle as a whole 
remains neutral, i.e., 2 = 0. When E p exceeds the field emission 
threshold, electrons are stripped (or field emitted) from the parti- 
cle, leaving the particle net positively charged, i.e., Q > 0. Similarly, 
when E p is below the field emission threshold, no electrons are emit- 
ted from the lower potential plate. However, when E p exceeds the 
field emission threshold, electrons from lower potential plate begin 
to tunnel across the vacuum gap and electric discharge sets in. In the 
figure, field emission electrons are indicated by e~ . 



position. Conductor plates with such material composition 
generated electric discharge at applied electric field strength 
of E p w 3850MV-m _1 . In a vacuum gap separated plane- 
parallel conductor plates, electric discharge is attributed to 
electrons from lower potential plate tunneling through the gap 
towards the plate with higher potential, thereby shorting the 
two plates and causing an electric discharge. In this regard, the 
phenomenon of electric discharge is intrinsically connected to 
the field emission threshold illustrated in Fig. Efc). Hence, 
in the aforementioned Zouache and Lefort's configuration, 
the field emission threshold of E cp < 3850MV • m _1 can be 
roughly approximated. This field emission threshold value 
may be compared with the electric field used by Szirmai to 
ionize his particle. In Szirmai's experiment, an electric field of 
E p « 4.67 MV • m _1 was sufficient to produce a net charge of 
~4x 10~ 15 C on a core-shell structured spherical aluminum 
particle of 3 |im in diameter. For the particle in Szirmai's ex- 
periment, the field emission threshold of E cs < 4.67 MV • m _1 
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can be assumed. Comparing the two values, one finds 



E cp _ 3850 MV -m- 1 
~E7 S ~ 4.67 MV-m" 1 



:824, 



which result suggests that for the configuration illustrated in 
Fig. [3] the particle can be sufficiently ionized when Zouache 
and Lefort's conductor plates are used. Further, by limiting 
applied electric field to a value much smaller than E cp and, 
yet, much larger than E cs , the particle can be ionized in the 
absence of field emission electrons originating from the plate 
held at lower potential. This prevents the ionized particle from 
neutralizing. For instance, the electric field value of E p m 
0.5E cp or E p ?s A\2E CS is much smaller than E cp , but it is still 
much larger than E cs . 

Aforementioned process of charging by field emission is 
schematically summarized in Fig. [5] To illustrate the mech- 
anism, particular electrons in the particle's shell and the core 
regions are labeled ej~ and , respectively. Similarly, partic- 
ular electrons in the lower conductor plate are labeled e7 and 

. It is understood that particle as a whole is initially electri- 
cally neutral, implying Q = 0. That cleared, an applied static 
electric field of magnitude E p = ||EJ| , where 



Ep ^> E cs , 



but E p <C E cp < E ds 



(14) 



is produced by connecting two plane-parallel conductor plates 
to a battery. Since E p ^> E cs , where E cs is the field emission 
threshold for the particle's conductor core (see Fig. |4j, 
tunnels through the dielectric shell and the particle as a whole 
becomes net positively charged. Because E p <C E dS) where 
E c / S is the field emission threshold for the particle's dielectric 
shell, no electrons are physically stripped from the shell and 
e\ remains confined to the shell, leaving the shell electrically 
neutral. Similarly, because E p <C E cp , where E cp is the field 
emission threshold for the lower conductor plate, no electrons 
can escape the surface of the lower conductor plate. The elec- 
trons, i.e., and , may redistribute themselves, but can- 
not physically escape the conductor plate's surface; which is 
the reason why cannot tunnel through the particle's dielec- 
tric shell to neutralize the positively charged core. To do so, 

must first escape the plate's surface, which is not possible 
since E p -C E cp . In conclusion, once the particle is charged 
and the electric field satisfies the condition defined in Eq. fl4l . 
the core-shell structured particle remains positively charged. 

Figure [5] presents a way to prepare yet another initial con- 
figuration, which is different from the one illustrated in Fig. [3] 
for the differential equation of Eq. < TT~0b . I shall briefly outline 
how such initial configuration might be prepared from Fig. [5] 
Starting with initially uncharged core-shell structured particle, 
the particle is ionized following the scheme illustrated in Fig. 
|5] Once the core-shell structured particle has been sufficiently 
ionized, the entire system is physically flipped over. Don't 
worry about the particle falling down because it won't. In 
fact, the positively charged particle sticks to the surface of the 
conductor plate held at lower of the two potentials. The rea- 
son for this is explained later in this section when the types of 
forces involved in the system are discussed. With the config- 
uration illustrated in Fig. [5]flipped over, the initial conditions 
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Figure 5: (Color online) Schematic illustration of ionization process 
of core-shell structured nanoparticle. For the illustration purpose, 
particular electrons in the particle's shell and the core regions are 
labeled and , respectively. Similarly, particular electrons in the 
lower conductor plate are labeled e7 and . Initially, the particle as 
a whole is neutral and Q = 0. After the field emission, the particle as 
a whole becomes net positively charged, Q > 0. 



are now specified by 



Zd 



(0)=0 and z d (0)=0. 



In the flipped over configuration, the applied electric field E p 
is in the e z direction. Now, to make this configuration suit- 
able for the differential equation of Eq. (flOt . the direction of 
applied electric field must be reversed. Thus, reversing the 
direction of E p to e z direction, this configuration becomes 
identical to the initial configuration illustrated in Fig. [3] ex- 
cept now Zd (0) = instead of Zd (0) = 0.25/z. Either choice is 
good for the initial condition of Zd(t) . That said, I shall keep 
using the initial conditions specified in Eq. ( fT3l and the initial 
configuration illustrated in Fig. [3]for the rest of this work. 

Returning to the plotting of Zd from Eq. ( fTOb with initial 
conditions specified in Eq. i\3i . the particle position as func- 
tion of time has been plotted in Fig. [6] using the parameter 
values defined in Eq. ( fTTb . The upper electrode is located 
at Zd = 0m in the plot. As it can be observed from the plot, 
the core-shell structured charged-particle executes an oscilla- 
tory motion between the two plane-parallel electrodes; and, 
such motion does not involve charge exchanges. One knows 
this because the governing equation of motion, i.e., Eq. ( fTol ). 
has been derived without any assumption of charge exchange. 
Such oscillatory behavior is fundamentally different from the 
traditional picture, which process assumes charge exchange 
mechanisms. The core-shell structured charged-particle re- 
bounds at Zd ~ 2.5 x 10~ 4 m from the upper conductor plate. 
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This rebounding position is too far from the lower conductor 
plate to account for any charge exchange processes even in the 
physical situations. In the plot of Fig. [6] the lower conductor 
plate is located at z c i ~h = 0.001 m. 
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Figure 6: Particle distance from the surface of upper conductor plate 
as function of time. For the plot, the values denned in Eq. d I 1 b have 
been used with initial conditions specified in Eq. d 1 3 b - Vt — Vl = 
8kV implies an applied electric field of E p = 8MV • m . The upper 
electrode is located at Zd = m. 

Hereafter, it is understood that charged-particle oscillation 
does not involve any charge exchanges for the rest of this pa- 
per. That said, what's the criterion for charged-particle os- 
cillation? Oscillatory modes occur when the effective charge 
carried by the positively charged particle satisfies the condi- 
tion given by 

where 



Z d,m (h - Zd,m) 

y(b 3 -a 3 )-b 3 y(b 3 -a 3 )-b 3 

w = ^-^ — + \ V~ 8<0 ' 

Z d,m [n - Zd,m ) 

and 

= 3k^ 

7 (K 2 + 2K 3 )bl+2(K 2 -K 3 )a3 < 

Similarly, for a negatively charged core-shell structured parti- 
cle, the oscillatory criterion is given by 

where 

1 1 

rj = j _ — > °- 

(h - Zd,m) Z d,m 



The mass dependence in Eqs. ([]]) and (O reveals that par- 
ticle with larger mass requires significantly larger effective 
charge compared to the particle with smaller mass to initiate 
oscillatory motion. Initially neutral particle can be charged or 
ionized by exposing it to a strong static electric field. For a 
nanoparticle, large portion of atoms composing it participate 
in the ionization, yielding in relatively large charge density 
per mass. However, for a macroscopic particle or an object, a 
great portion of atoms composing it does not participate in the 
ionization process and only those near the surface participate 
in the ionization due to electric field shielding effects. As a 
result, macroscopic particles have relatively small charge den- 
sity per mass. One may argue that the strength of electric field 
can always be increased to completely ionize the macroscopic 
particle. That, however, is not possible because, even in vac- 
uum, electric breakdown sets in at some point and everything 
neutralizes.— 

Based on this argument, the oscillation criterion specified 
in Eq. (Q]l is more likely to be satisfied by microscopic or 
smaller particles than by macroscopic counterparts. This im- 
plies the charged-particle oscillation presented in this paper 
is more likely to be observed from nanoparticle systems than 
from systems involving macroscopic particles. I shall now 
discuss why microscopic particles are more likely to satisfy 
the criterion of Eqs. (Q]) or (0 than the macroscopic counter- 
parts. To demonstrate this, I shall consider an aluminum ball 
of radius b = 1 .5 jim representing the smaller particle and an- 
other one with radius b = 3 mm representing the larger coun- 
terpart. To keep matters simple, I shall assume that the space 
between electrodes is a vacuum. 

In vacuum, jc 3 = 1 and KSqK 3 =2.78 x 10 _11 N _1 - rar 2 -C 2 , 
and the gravity constant is g = 9.8 m ■ s~ 2 . The mass densities 
of aluminum and aluminum oxide are p m \ = 2700 kg • m~ 3 
and p„, 2 = 3800kg m -3 , respectively. Thus, for a core radius 
of a — 1 .5 fim and shell thickness of b — a — 4nm, the total 
mass of the particle is obtained using Eq. (fT2b to yield m sa 
3.86 x 10~ 14 kg. The dielectric constant for the particle's shell 
is K"2 = 6. For the value of Zd.m, I shall choose Zdm — 0.25/z, 
where h = 1 mm. For the applied electric field, I shall choose 
E p = 8kV ■ m~' . Insertion of these values into Eq. (Q]) yields 

QT,mic = QT >5X10- 13 C, (15) 

where the notation Qr.mic denotes the microscopic particle. 

Now, I shall compute the same for the macroscopic coun- 
terpart. An aluminum ball of core radius a = 3 mm and shell 
thickness of b — a = 4nm has a mass of m w 3.05 x 10~ 4 kg, 
where Eq. (TT2l) has been used to compute the mass. To make 
sure this aluminum ball has sufficient room between the elec- 
trodes for oscillation, the gap between the two electrodes is 
increased to a value of h = 1 m. Keeping all other values same 
as previous, Eq. ([TJ gives 

Q T ,,nac = QT >9X10- 7 C, (16) 

where the notation Qr.mac denotes the macroscopic particle. 
Is this an experimentally obtainable value? The answer to this 
is maybe. It depends on what kind of electrodes are being 
used. Even in vacuum, one cannot increase the strength of 
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electric field indefinitely without electrical breakdown setting 
in, beyond which point everything neutralizes. 21 Comparing 
the two results, Eqs. ( TT3T > and ( [ToT l. Qr,mac is greater than 
Qr.mk by factor of a million. This result alone shows that 
the kind of charged-particle oscillation mechanism presented 
here, i.e., one that does not involve charge transfer processes, 
is most likely to be observed from microscopic or smaller par- 
ticles than from macroscopic counterparts. 

Because the system illustrated in Fig. |2ja) involves 
charged-particle executing an oscillatory motion, it radiates 
electromagnetic energy; and, the power of radiated energy can 
be obtained from Lienard radiation formula, 



Prad = 



4 



With the explicit expression for 'id inserted from Eq. (fTOb . this 
becomes 



Prod 



3c 3 

V 



Am 



{zc+bf 



[ 7 (b 3 -a 3 )-b 3 ]E p 



{h-zd-by 
[y(b 3 -a 3 )-b 3 ]E p 
(h-z d -b) 3 



[Zd + bf 



8E„ 



(17) 



The profile of Lienard radiation power corresponding to the 
oscillating core-shell structured charged-particle illustrated in 
Fig. [6] has been computed using Eq. ST% . The result shows 
train of emitted radiation power as the particle oscillates, as 
illustrated in Fig. [7] 

Radiation power (h = 1 mm, o, = 0.014 C/iri ) 




5e-06 le-05 

Time, t [si 



Figure 7: The profile of Lienard radiation power corresponding to 
the oscillating charged-particle illustrated in Fig. [6] 

It is worthwhile to discuss the types of forces involved in 
the plot of Fig. [6] The force responsible for generating such 
particle motion is Eq. ([8]l (or Eq. ©), of course. Since the 
force contribution from gravity is negligible in the oscillation 
regime, I shall only work with Fi and F2 of Eqs. © and ||7), 



respectively, for the discussion here. The force Fi is com- 
posed of the following three force contributions, 



Fl — f 1 1 +fl,2 + fl,3j 



where 



f 1 ,1 =e. 



1 Q T v 

16 s 2 ' 



. 1 Qt [y(ft 3 - a 3 )-ft 3 ] F 
fl ' 2 - e: l6 ? Ep ' 



fl,3 = -tz—Ep. 



From Eqs. (f3]l and it can be seen that 

V > and 1 > 7 > 0. 

This implies, 7 (b 3 — a 3 ) — b 3 < and the previous forces can 
be expressed as 



fl,2 
fl.3 ' 



1 



S 

-e 7 E 



Physically, f 1 1 represents the force between the charged par- 
ticle and the image charge formed on the surface of conduc- 
tor plate. Since the image (or induced) charge has opposite 
polarity, this force attracts the charged particle towards the 
plate. And, since fi 1 is in the direction of e z , this confirms 
such action. The f 1 3 physically represents the force by elec- 
tric field on the charged particle. Such force always pushes 
a positively charged particle in the direction of electric field. 
And, since f 1 3 is in the direction of — e z , which is the direc- 
tion of electric field, this also confirms such property. The 
remaining force term, f 1 2 , is a direct consequence of having 
a core-shell structured charged-particle. To be more accurate, 
f 1 .2 is a consequence of having particle with structure, which 
is not a "point" particle. Particles with structure can be po- 
larized by applied electric field and such property gives rise 
to fi 2- Consequently, this force vanishes in the absence of ap- 
plied electric field. For a positively charged particle, this force 
is induced in the same direction as the applied electric field. 

Similarly, the force F2 of Eq. (0 can be decomposed into 
the following three force contributions: 



F2 = f2,l +f2,2+f2,3j 



where 



12,1 = 



1 Qtv 



16 (h-s) 



2 ' 



1 Q T [y(b 3 -a 3 )-b 3 ] 
12,2 = e z — n 1 



'2,3 — 



16 

-ez—E p . 



{h-s) 3 



pi 
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Since y (fr 3 — a 3 ) — b 3 < 0, these can be expressed as 

f 1 
i2.i ~ -e 2 -5-, 

h.2 e z -r , 

Physically, f 2 1 is the force between the image charge and the 
charged particle. Since {2,1 is the force arising when charged 
particle is close to the lower conductor plate, this force must 
be directed towards the lower conductor plate. Since {2 1 is di- 
rected in — e z , this confirms such requirement. The {2,3 is the 
force on the positively charged particle due to the presence of 
electric field. Since electric field is in the — e z direction, so 
is f^3, as it must. Lastly, {22 is a consequence of having a 
charged-particle with structure and not a point particle. Par- 
ticles with structure can be polarized by applied electric field 
and such property gives rise to (22- Again, this force is al- 
ways in the direction of applied electric field for a positively 
charged particle. 

So, what is responsible for charged-particle oscillation? 
The second force, i.e., F2 of Eq. (0, cannot be responsible for 
charged-particle oscillation because {2 1, f2 2> and f2,3 are all 
directed in the — e z direction. However, the first force, i.e., Fj 
of Eq. ©, can give rise to charged-particle oscillatory modes. 
This is because force Fj contains which force points in 
the opposite direction of f 1 2 and fi 3. It is this competition 
between f 1 1 and the other forces, i.e., fi 2 and fi 3, that puts 
core-shell structured charged-particle in an oscillatory motion. 
Such mechanism is schematically illustrated in Fig. [8] In re- 
gion A, the dominant force is fi 2 and the magnitude of this 
force falls off like ~ 1 /s 3 with distance, s. However, in region 
B, contributions from f 1 2 weakens rapidly with distance and 
the force is dominated by f 1,1 , which force's magnitude falls 
off like ~ 1/s 2 with distance. Then, assuming the positively 
charged core-shell structured particle satisfies the oscillation 
criterion specified in Eq. ([T), the particle cannot escape re- 
gion B and enter region C, which region contains no oscilla- 
tion modes and any particles initially in such region ends up 
sticking to the lower conductor plate, as schematically illus- 
trated in Fig. [8] 

That said, a positively charged core-shell structured parti- 
cle initially in region B (or region A) would oscillate; and, 
the trace of such oscillatory motion over time would be rep- 
resented by the path 1 illustrated in Fig. [8] Here, the path 1 
represents the plot of Zd (t) versus time graph, where the time 
parameter is the horizontal axis. On the other hand, the same 
particle initially in region C would not have any oscillations, 
but it would follow the trace of the path 2, where the path 2 
represents the plot of Zd (f ) versus time graph with horizontal 
axis being the time axis. 

The sharp cusp in region A can be explained from the fact 
that fi 2 ~ — (l /s 3 ) e z is an extremely short range force. Since 
the magnitude of such force goes like ~ 1 /s 3 , it can generate 
very large impulse over short time. Since fj 2 is in — e z direc- 
tion, the particle is repulsed from the surface with very large 
force occurring over very short period. However, in region B, 



Mechanism for charged-particle oscillation 




lower conductor plate 



Figure 8: (Color online) Different forces dominate in each regions. 
In region A, the dominant force is fi 2- In region B, the dominant 
force is fj 1. In region C, the dominant force is F2. The path 1 and 
path 2 represent the plots of Zd (t) versus time graph, where the time 
parameter is the horizontal axis. Here, Vt >Vl. 

this force decays extremely rapidly and f^j ~ (1/s 2 ) e z domi- 
nates there. The direction of f 1 1 is in e z ; hence, the particle is 
pulled back to the upper conductor plate. This process repeats 
itself, resulting in an oscillatory motion. 

What happens when the magnitude of applied electric field, 
E p , is increased? The force f^i ~ (l/s 2 )e z is independent of 
E p ; therefore, lg does not change in Fig. [8] On the other hand, 
the force fi,2 ~ — (l /s 3 ) e z has an explicit dependence on the 
applied electric field, i.e., 

fi,2 ~ -e z ^-; 

therefore, the l& in Fig. [8] gets increased to l& + Al^ with 
increased E p . Here, I a and Zg represent the locations of bor- 
derlines for regions A and B, respectively. The result is that 
charged-particle trapped inside the region B is now forced 
to rebound more frequently due to the fact that the width 
of region B has been decreased by A/^. Consequently, the 
frequency of charged-particle oscillation increases with in- 
creased applied electric field; and, this is schematically illus- 
trated in Fig. [9] where it shows that new path 1 has higher fre- 
quency than the old path 1 . To validate this, Eq. ( fTOb has been 
plotted using the identical settings used to obtain the result in 
Fig. [6] This time, however, the strength of applied electric 
field has been increased from 8 MV -m _1 to 12MV m . The 
result is shown in Fig. [10] This result can be compared with 
the one in Fig. |6j which was obtained for E p = 8MV ■ m _1 . 
For the case of E„ = 12MV • m _1 , the frequency of charged- 
particle oscillation has been approximately doubled compared 
to that of the case of E p = 8MV • irr 1 . Notice that although 
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the oscillation frequency has been approximately doubled, its 
amplitude has been nearly halved. This must be so because 
the particle's oscillation frequency had been increased as a re- 
sult of reduced width of region B (or increased width of region 
A). 

Now, one cannot increase E p indefinitely to obtain higher 
oscillation frequencies because, eventually, the width of re- 
gion B would become zero. And, beyond that point, the 
charged-particle enters the region C and ends up sticking to 
the surface of the lower conductor plate. In region C, there 
are no oscillatory modes so any charged particle in that region 
gets attracted to the surface of lower conductor plate and stays 
there indefinitely. 

Mechanism for charged-particle oscillation 
(with increased applied electric field) 

upper conductor plate 



Particle position vs. time (h = 1 mm, CTj = 0.014 C/m ) 
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lower conductor plate 



Figure 9: (Color online) The width of region A has been increased 
to Ia + AIa as a result of increased strength of applied electric field, 



E„ 



—e z E p . The new path 1, old path 1, and path 2 represent the 



plots of Zd {t ) versus time graph, where the time parameter is the 
horizontal axis. Here, Vt > Vl- 

Having explained the kind of forces involved in the 
charged-particle oscillation, i.e., plots generated in Figs. © 
and ( fTOl ). it is now clear why there are cusps in the plot. The 
shape of these sharp turning points can be deceiving because 
these points are not really what they appear to be. In fact, 
these turning points are smoothly varying points and this is 
illustrated in Fig. [TT] where one of such sharp points has been 
enlarged for a view. At these points, the magnitude of force 
experienced by the particle falls off with distance like ~ 1/s 3 , 
where s is the distance between the particle's center of mass 
and the rebounding plate's surface. At very short separation 
distances, this repulsion force becomes extremely impulsive 
over very short period. But, nevertheless, ~ 1/s 3 is still a well 
behaved function because s cannot become zero, as the par- 
ticle cannot touch the surface of rebounding conductor plate. 
Doing so would require an infinite energy, which is not possi- 
ble. 
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Figure 10: Particle distance from the surface of upper conductor plate 
as function of time. For the plot, the values defined in Eq. i ll It have 
been used. Here, however, the value of Vj — Vi has been increased 
to Vj — Vi = 12kV, implying an applied electric field strength of 
E n = 12MV-m-'. The upper electrode is located at Zd = 0m. 



Particle position vs. time (/? = 1 mm, Oj = 0.014 C/rn ) 
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Figure 11: Particle distance from the surface of upper conduc- 
tor plate as function of time. The first sharp turning point near 
Zd ~ 0.000205 m in Fig. [T0]has been enlarged for a view, which 
shows a smoothly varying curve. 



What happens when the core-shell structured particle is 
negatively charged? In that case, the Fi and F2 of Eqs. (O 
and (0 get modified as 



Fi =e- 



127 



and 



F 2 = e, 



\Qt 



\v[ |y(fr 3 -a 3 )-fr 3 |£ p 
4s 2 4v 3 



\y(b 3 -a 3 ) -b 3 \E p 



4 (h - s) 



4(h-s 



2+ E P 



To distinguish the analysis here from the previous case involv- 
ing a positively charged core-shell structured particle, I shall 
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rewrite Fj and F2 as 
\Qt 



N, = e- 



|v| \y(b 3 -a 3 )-b 3 \E p 



As 2 



As 3 



and 



where 



N2 =n 2 ,i +n 2! 2 + n 2i 3, 



Igrj M 
16 (ft-j 



n 2 , 2 = e ; 



\Q T \ |y(fr 3 -a 3 )-fr : 

16 (h-sf " P ~^ Z (h — s) 3 ' 



-£"0 ~ e? - 



n 2 ,3 



12/ 



Since all of the terms are positive in force Ni , it cannot gener- 
ate any oscillations. On the other hand, the force N 2 contains 
both positive and negative terms; and, therefore, it can gener- 
ate oscillatory modes. Such is schematically illustrated in Fig. 

EI 

Mechanism for charged-particle oscillation 
(negatively charged-particle) 



upper conductor plate 




lower conductor plate 



Figure 12: (Color online) For the case of negatively charged-particle, 
the oscillation modes exist near the lower conductor plate, which is 
exactly opposite of the positively charged-particle case (see Fig. [8). 
In region A, the dominant force is n 2 . 2 ; and, in region B, the dominant 
force is 112,1- The path 1 and path 2 represent the plots of zd (f) versus 
time graph, where the time parameter is the horizontal axis. Here, 
V T > V L . 



except now <3\ = — 0.014C • m . Also, to account for the os- 
cillatory motion, the initial conditions have been specified as 



Zd 



(0)=0.75A and z,/(0)=0. 



Such initial conditions have been chosen because oscillatory 
modes only exist for Zd > 0.5h for negatively charged parti- 
cles. The result is plotted in Fig. Q~3] where it shows neg- 
atively charged core-shell structured particle oscillating near 
the lower conductor plate. In the plot, the lower conductor 
plate is located at Zd = 0.001 m and the upper conductor plate 
is located at Zd =0m. To show that charged-particle trajectory 
is represented by a well behaved function, one of the sharp 
turning points in Fig. [13]has been enlarged for inspection. As 
it can be seen in Fig. [T4J the cusp looking points are deceiv- 
ing because these are smoothly varying points. The magni- 
tude of force acting on the particle near these points falls off 
with distance like ~ 1/ {h — s) 3 . When the negatively charged 
core-shell structured particle is very close to the lower con- 
ductor plate, h — s becomes very small and this results in a 
very large force that acts to repulse the particle from the sur- 
face of the lower conductor plate. Nonetheless, this is a well 
defined force because h — s cannot become zero. Doing so 
would require an infinite energy, which is not possible. 

Particle position vs. time (h = 1 mm, o~ 1 = -0.014 C/rn ) 
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Figure 13: Particle distance from the surface of upper conductor plate 
as function of time. The lower electrode is located at id = 0.001 m 
and the upper electrode is located at Zd = 0m. The charged particle 
is negatively charged and it is oscillating near the lower electrode. 

This briefly summarizes the essence of this investigation. 
To complete the task, I shall now work out the detailed deriva- 
tions of key solutions used in this article. I shall begin by solv- 
ing the boundary value problem for the electrostatic potentials 
in regions Mi , M 2 , and M3 of Fig. [2] 



III. THEORY 



A. Free charge distribution 



To validate the argument illustrated in Fig. Q~2] Eq. (fTOb 
has been plotted using the same values specified in Eq. ( fTTT i. 



The correct specification of electric charge distribution is of 
crucial importance in any electrostatic boundary value prob- 
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Particle position vs. time {h = 1 mm, a 1 = -0.014 C/nT) 
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(a) Neutral spherical conductor in an otherwise 
uniform DC electric field E 



E 




(b) Charged spherical conductor in an otherwise 
uniform DC electric field E 



Figure 14: Particle distance from the surface of upper conductor plate 
as function of time. The first sharp turning point near id ~ 0.00086 m 
in Fig. [T3]has been enlarged for a view, which shows a smoothly 
varying curve. 



lem. When an uncharged, electrically neutral, spherical con- 
ductor is placed in an otherwise uniform electric field, the 
charges inside the conductor redistribute such that the poten- 
tial V\ is a constant there, as illustrated in Fig. [131 a). By def- 
inition, in an electrically neutral conductor, every charges are 
paired with one with opposite polarity. Therefore, the spher- 
ical conductor illustrated in Fig. [T"5l a). as a whole, is electri- 
cally neutral. 

Now, how do charges get distributed when you place a pos- 
itively charged spherical conductor in an otherwise uniform, 
constant, electric field? By definition, an electrically charged 
conductor has excess number of charges of one polarity that 
cannot be paired with one with opposite polarity. The paired 
ones do whatever they can to make net electric field zero in- 
side the conductor. The result is that the paired ones redis- 
tribute as illustrated in Fig. [T31a). What about the excess, un- 
paired, charges of same polarity? These must be redistributed 
such that the potential V\ is a constant inside the conductor. 
One such distribution, perhaps the only one, is illustrated in 
Fig- n"5l b). Assuming the charged spherical conductor has 
only the surface "free charges" and no volume "free charges," 
i.e., no excess charges embedded inside the volume, the free 
charges on the surface of spherical conductor must be uni- 
formly distributed over the entire spherical surface else the 
net electric field inside the spherical conductor would not be 
a zero. As an alternate explanation, the spherical conductor in 
Fig. [T31a), including its surface, represents an equipotential 
surface. When an excess of free charges of same polarity, say 
positive charges, is placed on such an equipotential surface, 
the charges get instantaneously redistributed over the surface 
due to Coulomb repulsion between the charges. The result is 
that these charges are uniformly distributed over the equipo- 
tential surface, as illustrated in Fig. [l31 b). 

That explained, I shall assume that surface "free charge" 
density <3\ in Fig. [2] is a constant which is uniformly dis- 
tributed over the surface r = a of the spherical conductor core 



E 




Figure 15: (Color online) Charge distributions of (a) electrically neu- 
tral spherical conductor and (b) a positively charged spherical con- 
ductor with surface free charge density, G\ , in an otherwise uniform 
electric field, E D . 



throughout this work. To generalize the problem, the dielec- 
tric shell surrounding the spherical core in Fig. [2] is allowed 
for a surface "free charge" density CTj at r = b. Introduction 
of free charge on the surface of dielectric shell is purely aca- 
demic. For realistic dielectrics, a% is negligibly small, if not 
zero. Hence, it can always be set to zero in the final solution. 
Therefore, I shall keep the problem simple by assuming that 
O2 is a constant which is also uniformly distributed over the 
surface r — b of the dielectric shell throughout this work. 

Lastly, although the illustration in Fig. Q3J considered only 
an excess positive charge case in which Qj > 0, the treatment 
throughout this work is not limited to such case only. The 
effective charge Qj can have either positive or negative polar- 
ities. 



B. Derivation of electrostatic potentials 

The apparatus for the problem is illustrated in Fig. EJa), 
where a core-shell structured charged-particle is placed be- 
tween two DC voltage biased plane-parallel conductors. Elec- 
trostatic potentials in regions M\ , M2, and M3 are described by 
Laplace equation, 



v 2 y 
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In spherical polar coordinate system, Fig. |2|b), Laplace equa- 
tion reads 



1 d ( 2 dV\ 1 d ( . dV 



= 0. 



r 2 dr \ dr J r 2 sin0 d6 \ dO 

1 d 2 V 



r 2 sin 2 9 d(j> 2 
For the system with azimuth symmetry, 
dV 



dtp 



o, 



the Laplace equation reduces to 

d ( 2 dV\ . 1 d ( . „dV 



sin0^— =0. (18) 



dr V dr J sin0 ^0 \ dd 
Equation ( fT8l has the general solution given by 

v{r,e) = i (V + ^V 

where coefficients A( and B/ : are constants, and Pg = Pf. (cos ) 
is the Legendre polynomial of order I. For regions M\, M 2 , 
and M3 in Fig. |2a), the electrostatic potentials are given by 



f=0 



(19) 
(20) 
(21) 



where coefficients A(, B(, Q, D^, and are to be joined to- 
gether by appropriate boundary conditions at the interfaces 
between regions. Equation ( fT9l does not contain terms like 
~ r~^~' because these terms blow up at the origin. 

Electrostatic potential inside of a conductor is constant; 
and, therefore, Eq. $1% becomes 

V l (r,9)=A = V , (22) 
where Vb is a constant. Potential is continuous at r = a, 

V 2 (a,d) = V l {a,d). 
From Eqs. ( f20b and ( f22b . it can be shown that 



or 



Bo - V + ^) P + £ [B e a e + ) P, = 0. 



(=1 



where Pq = 1. It follows that each Legendre polynomials (Po, 
Pi, P2, and so on) are linearly independent functions; and, 
therefore, the coefficient of each Legendre polynomials must 



be zero independently else this algebraic equation cannot be 
satisfied. Thus, 

B Q + Q )a - 1 =V , 
Beaf + CtcT 1 ' 1 =0. 
Solving for Co and Q yields 

Co = a(V -B ), 
Q = -B e a 2i+1 , l>\. 
From these results, Eq. ( f20b becomes 

aVo 



V 2 (r,0)=Bo(l-*) 

+ E B 



f=i 



r / r 

2i?+l 



(23) 



Equations (f2TT > and ( 1231 must be continuous at r = b, 

V 2 (b,6)=V 3 (b,e). 
With Eqs. ( l2TT i and ( |23l , it can be shown that 

JZi+l ■ 



fin 1 



aVp 
b 



E^K 



b t+i 



Do + ^ + iUb*' E < 



b t+\ 



or 



Bo 1-7 



— - D °-T 



E 



Be lb 



.,21+1 



b t+l 



DiV 



Po 



b e+i 



where Pq = 1 . Again, Legendre polynomials (Po, P\,Pi, and so 
on) are linearly independent functions; and, therefore, the co- 
efficient of each Legendre polynomials must vanish indepen- 
dently else this algebraic equation cannot be satisfied. Hence, 



En (a 
D + ^=B (l-- 



flVb 



D e b l + 



E, 



B ( b 



,2i+X 



b e+i \ " h e+i 
and the following coefficients are obtained: 

£0 = Bo (b — a) + oVq - D b, 
E ( =Bdb 2e+1 -a 2M )~D e b 2M 



£>l. 

Using these results, Eq. (|2TT> becomes 

b\ B (b-a)+aVo 



V 3 (r,0)=D 1 



E 



De r 



~J+r 



Be (b 2e + 



1 JU+Y\ 



r e+i 



(24) 
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Equation d24l i must simultaneously satisfy the boundary con- 
ditions at the surfaces of the upper and the lower conductor 
plates illustrated in Fig. Eta). 

In Cartesian coordinates, the surface of the upper conductor 
plate is described by the z = s plane and the surface of the 
lower conductor plate is described by the z = s — h plane. At 
distances sufficiently far from the particle (or very close to the 
surface of conductor plates), the potential inside of the parallel 
plates can be approximated as 



V p = - [ E p -e z dz' + V L , 

Js-h 



where Vl is the voltage applied to the lower conductor plate, e z 
is the versor along the Cartesian z axis, and E p is the electric 
field inside of the parallel plates in the absence of the charged- 
particle. The expression for E p is given by 



E, : 



-e*£(Vr-Vi), 



from which the V p can be obtained: 



V p (z)=E p (z-s + h) + V L , 



where 



E„ 



(|Vr-Vi|). 



(25) 



(26) 



(27) 



In spherical polar coordinate system, the Cartesian coordinate 
z is represented by 

z = rcos0 

and Eq. d26b becomes 

V p (r, 9) = E p (rcos 9-s + h) + V L . (28) 



For r very large, but not infinite in extent, the contributions 
and ~ r~'~' become negligible in Eq. 



from terms like 



(l24l and the V3 takes the form given by 

V 3 (r,9)^D +'^D i /p e , 



(29) 



where b -C r < °°. At distances sufficiently far from the parti- 
cle, V 3 (r,0) « V p (r, 9) . 

One may recall a typical problem in electrodynamics, 
wherein a charged sphere is immersed in an otherwise con- 
stant and uniform electric field. For instance, assuming uni- 
form electric field is along the e z axis, a useful boundary con- 
dition is that at infinity, electric field is just E = e z ||E|| rcos 9. 
Now, at distances which are infinitesimally close to the con- 
ductor plate's surface, electric field must be perpendicular to 
the plate's surface. This is because the surface of conductor 
plate is an equipotential surface and electric fields are perpen- 
dicular to the equipotential surface by definition, of course. In 
this regard, Eq. d28l mimics the electric field boundary con- 
dition at infinity for the textbook problem in electrodynam- 
ics, wherein a charged sphere is immersed in an otherwise 



constant and uniform electric field and the student is asked to 
solve for the potential around the sphere. 

That said, Eqs. d28l and d29b are equated to yield 



D + D,rcos9 



i=2 



& E p r cos 9 +E p (h-s) + V L . 

Matching the coefficients of the like Legendre polynomials 
yield 



Do 



;E p (h 
'-'pi 

i0, t>2. 



s) + V L , 



Using these results, Eq. ( l24b becomes 



V 3 (r,9)^E p (h-s) + Vi- 



E p r cos 9 



[B (b-a) + aV - bE p (h - s) - bV L ] - 



\B\ (b 3 -ai)-b 3 E p ] 



cos 9 



E 

e=2 



Be (b 



2t+\ 



(30) 



where it is understood that b <C r < °°. The electrostatic po- 
tential, which satisfies the Laplace equation, is a second order 
differential equation. Therefore, its derivatives must be satis- 
fied at the boundaries. The remaining unknowns, Bo, B\, Bg 
for I > 2, and Vq are evaluated from the statement about the 
discontinuity of electric displacement at r = b and at r = a. 

At r = b, the normal component of the electric displacement 
suffers a discontinuity given by 



[e r -D 3 (r,0)-e r -D 2 (r,0)] \ r=b = a 2 , 



(31) 



where 02 is the surface free-charge density at r = b, the e r is 
a unit vector pointing in the radially outward direction, and 
D2 and D3 represent electric displacements in regions M2 and 
M3, respectively. In the linear dielectric approximation, the 
electric displacement can be expressed as 



Vi(r,e) = -eoKiWi(r,0), 



(32) 



where fQ is the dielectric constant in region M, and £0 is the 
electric permittivity of the free space. Hence, Eq. OTT ) can be 
expressed as 



[jc 2 e r • VV 2 (r, 6) - K- 3 e r • VV 3 (r, 9)} \ r=b 



02 
£o' 



(33) 



In spherical polar coordinate system, the V operator is defined 
by 



d 



1 d 
* e -rd9 



1 



rsin 9 d(j) 



and Eq. d33l becomes 



dV 2 {r,9) dV 3 {r,9) 
K 2 ^ - K 3 ^ - 



dr 



02 



~i (34) 
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which constitutes the Neumann boundary condition at r = b. 
In explicit forms, the derivatives in Eq. ( 134b are evaluated as 



dr 



(£+l)a 



2t+l 



rt+2 



(35) 



and 



dV 3 (r, 0) 



E p cos 



[B (b-a) + aV Q - bE p (h - s) - bV L ] -j 

r L 



■2[B t (b 3 -a 3 ) -b 3 E p ] 



cos 



" (£+l)B e {b 2M -a 2M ) n 
■ L jTi Pl 

1=2 r 



(36) 



where Eqs. ( 1231 and (T30b have been used. Insertion of Eqs. 
d35l and © into Eq. OH yields 

— = {B [a (K 2 - K 3 )+bK 3 ] 

~a(K 2 ~ K 3 )V - bKiEp (h -s)- bK 3 V L } 
2a 3 * 



"2 1 l+-p- )filCOSf) 



K 3 ^[ Bl (b 3 -a 3 )- i 3 ^] - £ P | cos 



£ + 1 



yi+2 



k 2 £ 

b i-e 

2M (K 2 -K 3 )+b 2e+1 K 3 



or 



-a (k 2 - K 3 ) V - bx 3 E p {h-s)\ Pq 

+ ^2(1 + ^)^1 

+ K 3 1 1 [Bi (Z> 3 - a 3 ) - - £ p | Pi 



1=2 
£+1 

a 



k 2 £ 



b t+2 



b l 

2t+l 



(K 2 -K 3 ) + b 2i+l K 3 



where Pq = 1 and P\ = cos 9 . Because each Legendre polyno- 
mials of order £ are linearly independent functions, this alge- 
braic relation can be satisfied if and only if the coefficients of 
each Legendre polynomials vanish independently. Hence, 

B [a(K 2 -K 3 )+bK 3 ]-a(K 2 -K 3 )V 

-bK 3 [E p {h-s)+ V L ] - — = 0, 
£0 



1 [d 3 - a 3 ) -b 3 E p ] -E p 
3 



+ K - 1 l + ) Bl =0 - 



and 



(£+1 

B iw2 



a 2e+l {K 2 -K 3 )+b zt+L K 3 



K 2 £ 

ui-e 



= 0. 



One reads off immediately that 

bK 3 [E p (h - s) + V L ] + b 2 e 1 a 2 



Bo 



a(K 2 - K 3 )+bx 3 
a(K 2 -K 3 )V 



a(K 2 -K 3 )+bK 3 '' 



B\ = 



3K 3 b 3 E p 



(K 2 + 2K 3 )b 3 +2(K 2 -K 3 )a 3 '' 



and 



B £ = for £ > 2. 



(37) 



(38) 



(39) 



With coefficients Bo,B\, and B(> 2 defined, Eqs. ( 123b and (f30b 
become 



V 2 (r,e)nB 1 



and 



r 

+ Bi ^1 rcos0 

V 3 (r, 0) « E p (h -s) + V L + E p rcos 9 

+ [B Q (b-a) + aV - bE p [h - s) - bV L ] - 



(40) 



[ Bl (b 



3 a 3 ) 



-b 3 E p ] 



cos 



(41) 



where Vo is the only unknown. 

The Vo is evaluated from the statement about the disconti- 
nuity of electric displacement at r = a. At r = a, the normal 
component of the electric displacement suffers a discontinuity 
given by 



[e,--D 2 (r,0)~e,--D 1 (r,0) 



01, 



(42) 



where (J\ is the surface free-charge density at r = a and Di 
is the electric displacement in region M\ . Repeating the same 
procedure outlined from Eq. OTb through Eq. d34l l. it can be 
shown that 



„ dVl(r,0) dV 2 (r,0) 

K\ - K 2 - 



dr dr 
Since region M\ is a conductor, 
dVx (r, 9 



£o 



= 
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and Neumann boundary condition at r = a becomes 



dV 2 (r,0) 



£0*2' 



(43) 



Using the results in Eqs. (f37l >. (l38~l l. and ( |39l l, the derivative in 
Eq. d43l is readily computed from Eq. d35l >. 



dV 2 (r,9) 



{Bo - V ) ■ 



-By 1 



2a 3 



cos 0. 



With this result, Eq. d43l becomes 
h3Bicos0 



(44) 



The cos 9 in Eq. d44l i can be eliminated by integrating both 
sides over the spherical surface at r = a, 

/Bq-Vq + ^ \ 
=o J>p=o \ a J 

-2k 

-a" sin9d9d<p, 



9=oJ(j>=0 e K 2 



yielding 



Bo-V = - 



BoKz' 



(45) 



What I have just done here only surmounts to the computing 
of total free charge on the sphere of radius r — a. For instance, 
in Eq. (1441) . one can integrate both sides over the surface r = a 
of a sphere. The right hand side yields total free charge on the 
surface r = a, ignoring the extra constant factor. The left hand 
side yields terms with cos 9 eliminated, as this term has been 
integrated over. Canceling out the common terms yields Eq. 

<E3. 

That explained, Bq is inserted from Eq. fl37T i into Eq. (|45]> 
to solve for Vb; and, this yields 



Vq = V l - 



' — a)a\ 



a 2 <3\ 



■b 2 a 2 



bSoK 2 beoK 3 
+ E p (h-s). (46) 

With Eq. (|46T >. the coefficient Bq of Eq. (|37| | becomes 



Bq = Vl 



a(2b — a)G\ a 2 Gi+b 2 (J 2 



be K 2 
+ E P (h-s). 



bE K?, 



(47) 



With coefficients B\, Vo, and Bq defined respectively in Eqs. 
d38l l. d46l i. and (|47] >. the electrostatic potentials for regions Mi , 
M 2 , and M3 are obtained from Eqs. d22b . d40l i. and ( fTTb . They 
are 



Vi=Vi + a+E p (h-s), r<a, 



(48) 



V 2 (r,d) = V L + P+Ep(h-s + yrcosd) 
X a 3 jE p cos 9 
r r L 



a < r < b, 



(49) 



and 



V 3 (r,6) =V L + E p (h-s + rcos9) + - 



[y(b 3 ~a 3 ) -b 3 ]E p cos G 



-C, r>b, 



(50) 



where a, /J , 7, k , and v are defined as 

a{b — a)(5i a 2 0\+b 2 2 
be K 2 be Q K 3 ' 

a(2b — a)0\ a 2 0\+b 2 o 2 



a 
P 

r= 
x = 

V = 



be K 2 be K 3 

3K 3 b 3 

(K 2 + 2K 3 )b 3 + 2(K 2 -K 3 )a 3 ' 

a 2 0\ 



(51) 



£0^2 

2a{b — a)<j\ a 2 0\+b 2 2 



60*2 



eoJc 3 



For all of the treatment hereafter, only the derivatives of V 3 , 
in particular, the normal derivatives associated with the plane- 
parallel plate electrodes, are of importance. Therefore, the 
explicit expression for the constant C in Eq. < f30b is not of 
much concern here. 



C. Induced surface charges on conductor plates 

In spherical polar coordinate system, V operator is defined 

d 



by 



d Id 

V = e ^ +ee 7^ 



1 



rsin 9 d(j) 



where 



e r = e T sin 9 cos <p + e y sin 9 sin <j) + e z cos 9 , 
eg = e x cos 9 cos + e y cos 9 sin <p — e z sin 9 , 
= — e x sin <j> + e y cos . 

Hence, the e z component of V operator is given by 

. d sin 9 d 

e, (e, ■ V) = e,cos0- — e, — . 

dr r 69 

Using the form defined in Eq. d32i l. the electric displacement 
in region M 3 is given by 

D 3 (r,e) = -eoK 3 VV 3 {r,0). 

The e z component of D3 (r, 9) is obtained by replacing the V 
with the e z (e z ■ V) operator and this gives 



D 3;z (r, d) = EoK 3 e z 



sin9dV 3 (r,9) 
~r d~9 



cos 9 



dV 3 (r,9) 



where the notation D 3 - Z (r,9) denotes the e, component of 
D 3 (r, 6) . With V 3 (r, 9) of Eq. (O, the e, component of elec- 
tric displacement in region M 3 is given by 
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D 3;z (r,0) = £ fC 3 e z jocose 

+ \J\ J LJi( 3cos 20_ 1 



(52) 



The surface of the upper conductor plate is described by 
the Cartesian z = s plane. In the spherical polar coordinate 
system, the surface of the upper conductor plate is described 
by 

s 

cos = 



yx 2 + y 2 + s 2 

Insertion of the expression for cos 6 into Eq. d52l yields 

( 3[y(b 3 -a 3 ) -b 3 ] E p s 2 



(x 2 + y 2 +s 2 ) 5/2 



Vs- [y(b 3 -a 3 ) -b 3 ]E p 
(x 2 +y 2 +s 2 f 2 



(53) 



At the surface of the upper conductor plate, the electric dis- 
placement suffers a discontinuity given by 



(54) 



where Gi up is the induced surface charge density on the surface 
of the upper conductor plate and D MCp;z is the e z component of 
the electric displacement inside of the upper conductor plate. 
Since the electric displacement inside of the upper conductor 
plate is zero, Eq. ( T54T i reduces to 

e z -D ix (x,y,s) = -G iup 

and the surface charge density is given by 

r '3[y(b 3 -a 3 )-b 3 ]E p s 2 



Oiup = -£0*3 



(x 2 +y 2 + s 2 ) 5/2 



Vs- [y{b 3 -a 3 ) -b 3 ]E p 



(55) 



(x 2 +y 2 + s 2 ) 3/2 

where Eq. d53l has been inserted for D 3;z (x,y,s) . 

The surface of the lower conductor plate is described by the 
Cartesian z = s — h plane. In the spherical polar coordinate 
system, the surface of the lower conductor plate is given by 

s — h 

cos 9 



\Jx 2 +y 2 + {s-hf 



and Eq. d52b becomes 



e z £oK3 \ 



3 [y{b 3 -a 3 ) -b 3 ]E p (s-h f 



+ 



x 2 +y 2 + (s-h) 2 
V(s-h)- [y(b 3 -a 3 )-b 3 ]E p 



x 2 +y 2 + (s-hf 



3/2 



5/2 



-E„ 



(56) 



At the surface of the lower conductor plate, the electric dis- 
placement suffers a discontinuity given by 

e z • D 3;z {x,y,s- h) - e z • D /c . p;z (x,y,s - h) = a ap , (57) 

where Ou p is the induced surface charge density on the surface 
of the lower conductor plate and D/ cp;z is the e z component of 
the electric displacement inside of the lower conductor plate. 
Since the electric displacement inside of the lower conductor 
plate is zero, Eq. ( T57l i reduces to 

e z' D 3; z (^3','S-/l) = Gilp 

and, with Eq. ( T56l l inserted for D 3;z (x,y,s — h) , the induced 
surface charge density is given by 



Gii p = e fc 3 



3 [y(b 3 -a 3 )-b 3 ]E p {h-s) 2 



x 2 +y 2 +(h-sf 



5/2 



v(h-s)+ \y(b 



3 -a 3 ) 



-b 3 ]E p 



x 2 +y 2 + (h-sf 



3/2 



-E 



p ( » 



(58) 



where (s — h) has been re-expressed as — (h — s) purely for 
convenience. 

In the limit the parallel plates become infinite in extent, the 
total of induced charges on the surfaces of each conductor 
plates must add up to the total charge carried by the particle. 
To check on this, Eqs. (l55l l and (l58l are integrated over the 
surfaces of infinite parallel conductor plates with gap h. For 
convenience, I shall perform the integral in the polar coordi- 
nate system. In terms of the polar coordinates, Eqs. ( T55l l and 
become 



a iup (p,s) = -£oK 3 



'3[y{b 



3 -a 3 ) 



■b 3 ]E p s 2 



v (p 2 + s 2 f 2 
vs - [y(b 3 -a 3 ) -b 3 ]E p 
(p 2 + s 2 ) 3/2 



(59) 



and 



Gn p (p,s) = eqk 3 < 



3 [y{b 3 -a 3 )-b 3 }E p {h-s) 



p 2 + (h-s) 



5/2 



v(h-s)+ [y{b 3 -a 3 )~b 3 ]Ep 



p 2 + (h-sf 



3/2 



-E 



p ( > 



(60) 



where p = y/ 'x 2 +y 2 . Since the surface in polar coordinate 
system is symmetric about its axis, the total induced charges 
on both conductors can be performed as follow: 



QiT = Qiup + Qilp 
r-2n 



=0Jp=0 



{p,s) + Gii p (p,s)] pdpdtj) 



2k [ [G hip {p,s) + Gap (p,s)] pJp, 

Jp=0 
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where Qi up and Qn p are respectively the total induced charge 
on the surface of the upper and the lower conductor plates. 
With Eqs. (|59]l and d60j, the Q iT becomes 



Qt 



3[y(b 3 -a 3 ) -b 3 ]E p s 2 



JIt\,K, ./,. | (p 2 + S 2 ) 5/2 

3 [y(b 3 -a 3 ) -b 3 ]E p (h-s) 2 



p 2 + (h-s) 



5/2 



+ 



VS- [y^-a^-b^E,, 
{p 2 + s 2 f 12 

v (h - s) + [y(b 3 - a 3 ) - b 3 ] E p 

3/2 



p 2 + {h~sf 
Equation doTt involves the following integral types: 



pdp. (61) 



pdp 



(p 2 + C 2 ) 3/2 sfp 2 ^ 



(62) 



pdp 



1 



2 1 ^2^5/2 



(p 2 + c 2 ) 



3(p 2 + c 2 f 2 



1 

3~ 3 



TT- (63) 



With the integral formulas of Eqs. ( f62b and d63l ), the gir of 
Eq. doTt is integrated to yield 



Q/r = -47reoK:3V. 
Insertion of the explicit expression for v from Eq. ( BTT l yields 

Qff = -(G6+Gi + G2), (64) 



where 



g& = 87ra (b-a)0\ 

Qi =Ana 2 G\, 
Q 2 = \%b 2 a 2 . 



*3 

K 2 



The three quantities are identified as follow. The Q\ and Q 2 
are the "free charges" on the surfaces at r = a and r — b, re- 
spectively. The Qt, is the charge contribution arising from the 
presence of a dielectric shell surrounding the metallic core. 
This contribution vanishes in the absence of free charge on 
metallic core (i.e., Ci = 0) or dielectric shell (i.e., b — a = 0). 



D. Particle dynamics 

Two major electrostatic forces are acting on the core-shell 
structured charged-particle in Fig. [2] One such force is the 
electrostatic force between the induced charges on the surface 
of the upper electrode plate and the charged-particle. This 
force is denoted by Fi. The other force is the electrostatic 



force arising between the induced charges on the surface of the 
lower electrode plate and the charged-particle and this force is 
denoted as F2. The net force exerted on the charged-particle 
by induced charges on each surfaces of the conductor plates is 
therefore given by 



F = Fi+F 2 



jdE, + JdE 2 I , 



(65) 



where Qt is the effective charge carried by the charged- 
particle and dE\ and dE 2 are respective differential electric 
fields corresponding to the upper and lower electrode plate 
surfaces Si and Si, respectively. For instance, dE\ is the 
differential electric field associated with the induced surface 
charge at location Ri in Fig. |2fa). Similarly, dE 2 is the differ- 
ential electric field associated with the induced surface charge 
at location R2 of Fig. |2ja). 

The presence of extra factor of 1/2, the negative sign, and 
the exact form of Qt in Eq. d65l ) can be explained as follow. 
The extra factor of 1/2 in Eq. (|65] l comes from the fact that 
each parallel conductor plates sees only an hemisphere of the 
charged-particle. The effective charge carried by the particle 
is identical in magnitude to the giT of Eq. d64l ), but with op- 
posite charge polarity. Thus, 



Qi 



Q 



IT 



or 



Qt = Sna (b-a)o 1 — +4n (a 2 a x +b 2 o 2 ) 

K 2 y ' 



(66) 



The negative sign in Eq. (1651 1 is necessary for specifying cor- 
rectly the direction of the forces exerted on the core-shell 
structured charge-particle by induced surface charges from 
each parallel conductor plates. To demonstrate this, the in- 
tegrals in Eq. d65l l can be represented by 



dEi 



47C£3 



Ik 



P gRiPidpidfr 



oJ P ,=o (R, -R, 



>3/2 



(67) 



where Q is the induced surface charge at location R, (i.e., 
i= 1,2) in Fig. |2|a). Now, suppose if £1 is positive, then 
the direction of dE\ must be in — Ri, as it can be inspected 
from Fig. |2|a). On the other hand, if q\ is negative, then the 
direction of dE\ must be in Rj. The same argument can be 
said for those involving q 2l dE 2 , and K 2 . And, this explains 
the presence of negative sign in Eq. d65l ). 

That said, using the form defined in Eq. d67l ). the force 
expression of Eq. ( f65T > becomes 



F; 



Qt 



In 



q^iPidpidfyj 



87Ce 3 J^ = oJp i= o (R, R, 



3/2 



(68) 



where 2 = (1 ,2) , £i = <J iup of Eq. d55), q 2 = Gu p of Eq. d60l l, 
and £3 is the electric permittivity of the region M3 . The ex- 
plicit expression for R,-, which defines the position of the g 
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associated with dSi as illustrated in Fig. [2 a) for i = (1,2) , 
are given by 

Ri = e x pi cos <j>\ + e v pi sin ty\ + e z s, (69) 
R2 = e. r p 2 cos <p 2 + e v p2 sin tj) 2 + e z (s - h) , (70) 

where h > s. 

The force exerted on the particle by the induced charge on 
the surface of the upper conductor plate is obtained by insert- 
ing R[ of Eq. $69$ into Eq. fl68). This yields 



' 1 = 



+e v 



2 k 



giPicos^i 



giPisin^i gis 



2^3/2 



pidpid^i. (71) 



The two terms in the integrand with cos 0i and sin 0i vanish 
when integrated over d(j). Thus, Eq. $7T\ reduces to 



Fi = -e z - 



Q T s f p a iup pidpi 



4e 3 Jo ( p 2 + 5 



2 4.^V2' 



(72) 



where 51 = a, Mp . Insertion of the explicit expression for (7„ (/ , 
from Eq. d59l into Eq. ( |72] | yields 

Gr* f p \l>[y{b 3 -a 3 )-b 3 }E p s 2 

* 1 — — , — 



I (p 

vs- [y(fo 3 -a 3 ) -£ 3 ]£ P 



2 1 c 2^ 



(Pl+siy 



2^^3/2 



P\dp\- 



[Pt + s 



Equation ( 1731 involves the following type of integrals: 



f p Pidpi 


1 


1 


Jo (p2 + ,2) 4 


~ 6? 


6(p 2 + s 2 ) 3 


f p P\dp\ 


1 


1 


Jo (p2 + ,2)3 


~ 4s 1 


4(p 2 + *2)2 


f p P\dpi 


1 


1 


'0 (P? + ^ 2 ) 3/2 


5 


\/p 2 + s 2 



(73) 

(74) 
(75) 
(76) 



Insertion of Eqs. <l74l . ( 1751 . and d76j into Eq. (|73]l yields 

2r / V Vs 2 [y{b 3 -a 3 )^b 3 ]E p 



Fi =e- 



16 1 S 2 (p2 + iS 2)2 5 3 



[ 7 (* 3 



-fo 3 ]£ p 5 2[y(fo 3 -a 3 ) -fc 3 ]£ p ,s 3 



AEpS 



v / p 2 + s 2 



(p 2 + s 2 Y 
4E n 



(p 2 + s 2 ) 



(77) 



where b <s <h — b. Equation (fTTJ is the force exerted on the 
charged-particle by the induced charges on the surface of the 
upper conductor plate. 



The expression for the force exerted on the particle by the 
induced charges on the surface of the lower conductor plate is 
obtained by inserting R2 of Eq. d70l into Eq. d68l l. Repeat- 
ing the similar procedure outlined in Eqs. (|7T1 and d72l . one 
obtains 

Q T {h-s) f p Ouppidp! 



4e 3 



p 2 + (h- s y 



3/2 ■ 



(78) 



Insertion of the explicit expression for Gn p from Eq. (l60l > into 
Eq. (ED yields 



F 2 =e : 



Qt (h 



3 [y(b 3 ~a 3 ) -b 3 ]E p (h-s) 2 



pi + (h-sy 



V(h-s) + [y(b 3 -a 3 ) -b 3 ]E p 
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Using the integral formulas from Eqs. ( 1741 . ( 1731 ). and ( 1761 ) 
with s replaced by h — s, Eq. ( 1791 becomes 



F2 = e *16< 
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where b <s < h — b. Equation (f80b is the force exerted on the 
charged-particle by the induced charges on the surface of the 
lower conductor plate. 

For a parallel plate system, which is microscopically large, 
but macroscopically small, the forces in Eqs. ( 1771 1 and ( f80b 
can be approximated by making p go to infinity. This ap- 
proximation is certainly valid for very small charged-particles 
confined between large parallel conductor plates. In the limit 
p goes to infinity, Eqs. (fTTJ and (fSOb simplify in form as 
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where b < s < h — b and v is defined in Eq. ( BTT i. 

1a{b — d]0\ a 2 (J\+b 2 <77 

v = — 1 . 

e K- 2 £oK- 3 

Notice that the resulting forces in Eqs. dSTI ) and d82l > are 
now just one dimensional forces; that is, Fi = F[ (s) and 
F2 = F2 (s) , where the parameter s is the relative distance 
between the center of mass point of core-shell structured 
charged-particle and the surface of upper conductor plate. The 
dynamics of charged-particle system illustrated in Fig. |2]has 
now reduced down to solving a nonlinear ordinary differential 
equation. 

Insertion of Eqs. (T8TT > and (l82b into Eq. d65l > yields the total 
force exerted on the charged-particle by the induced charges 
on the surfaces of parallel plate conductors. The result is 
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If the gravitational effect is included, the force experienced by 
the particle is 



or 
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(83) 



where m is the mass of the particle, g = 9.8 m • s~ 2 is the grav- 
ity constant, and the gravitational force has been assumed to 
be in the — e z direction. Since Qt is related to V by 

Q T =47z;£ok-3V, (84) 

the Ft of Eq. ( f83l may be re-expressed for convenience as 
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The dynamics of oscillating charged-particle is given by 




(86) 



where c = 3xl0 8 m-s 1 is the speed of light in vacuum. The 
left hand side of Eq. d86l l can be differentiated to give 




where v = dv/dt. Knowing that 
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equation (|87| > becomes 
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Multiplying both sides by c 2 (l — v 2 /c 2 ) 3 ^ 2 yields 

m \ c L ) 
Since v = s and v = s, Eq. d88l becomes 

se 7 = — 1 — =• 



(88) 



where s = d 2 s/dt 2 . With F^ explicitly inserted from Eq. 
the expression for s becomes 
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where e z has been dropped for convenience. It is convenient 
to re-express Eq. d89b in terms of the variable z c i illustrated in 
Fig. |2ja). Two variables, s and Zd, are related to each other by 

s = Zd+b, s=id, s = Zd, (90) 
where b is a constant. Hence, in terms of Zdi Eq. d89l becomes 
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(91) 



Equation j9TT l governs the dynamics of an oscillating charged- 
particle, subjected to high electrostatic fields, at all speeds. 



E. Criterion for charged-particle oscillation in the absence of 
charge transfer process 

The criterion for charged-particle oscillation in the absence 
of charge transfer process between the rebounding electrode 
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and the charged-particle can be obtained by analyzing the 
force expression of Eq. ( |85l >. 



F T = e. 




(92) 



The kinematics of charged-particle motion associated with the 
force of Eq. d92| i is illustrated in Fig. [6] where the param- 
eters s and Zd are related by s = z c i + b. For the plot illustrated 
in Fig. [6] the turning points of charged-particle motion occurs 
approximately at z c i = Zd.m — 0.25/z. 

By definition, when the particle is in vicinity of the turn- 
ing point, but not past it, the net force acting on the particle 
is directed in the opposite direction of particle's motion. Im- 
mediately past the turning point, the particle's motion is in the 
same direction as the net force. Thus, immediately past the 
turning point, the net force satisfies 

Fr (zd, n ) > 
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Referring to the plot in Fig. |6l this result simply states that 
net force is in e z direction at Zd = Zdm- This relation can be 
rewritten as 
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Notice that < 7 < 1; therefore, y/ < 0. Thus, the previous 
relation can be rewritten as 
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Solving for E p , I obtain 
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In terms of the effective charge Qj defined in Eq. d84t . this 
result becomes 
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(94) 



By definition, E p > because it is the magnitude of elec- 
tric field. Therefore, the right hand side of Eq. d94l i must 
be positive. But, according to the plot illustrated in Fig. [6] 
Zd.m < b/2; and, thus 



> 



- Zd., 



For a positive particle, Qj > and Eq. ( f94b can always be sat- 
isfied for a sufficiently ionized positive core-shell structured 
particle. For convenience, I shall express Qj in terms of E p . 
Equation ( |94| l can be rearranged to yield 
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Utilizing the quadratic formula, this can be solved for Qj to 
yield 
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When the core-shell structured charged-particle is nega- 
tively charged, how does the charged-particle oscillation crite- 
rion get modified? For the negatively charged core-shell struc- 
tured particle, v < 0, where 
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Thus, the Ft of Eq. (|92l ) becomes 
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Furthermore, since y(b 3 — a 3 ) — /? 3 < 0, 1 have 
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Now, immediately past the turning point at Zdm = 7.5 x 
10 _4 m in Fig. [13] the net force acting on the particle satis- 
fies the condition given by 
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where Eq. d97l ) has been substituted in for the explicit expres- 
sion. This relation is rearranged as 
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where Eq. d84l has been used. Thus, E p becomes 

,2 



Amg 



E„< 



1 



1 



1 



(h-Zd,m) 2 z d,m 



\6mg 



(98) 



Now, E p cannot be negative because it is the magnitude of 
electric field. Since Zdm > h/2 in Fig. [13] it must be true that 
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Therefore, Eq. d98l ) can always be satisfied for a core-shell 
structured particle which is sufficiently charged negatively. 
Rearranging Eq. d98l ), I have 
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Introducing r\ , 
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Utilizing the quadratic formula, I find 
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Since \Qj\ cannot be negative, one must make sure the right 
side is positive. But, since 
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Equation d99l is the oscillation criterion for the negatively 
charged core-shell structured particle. The resulting form is 
different from the positively charged case, i.e., Eq. ( |95T l, due 
to the fact that negatively charged particle oscillates near the 
lower conductor plate whereas the positively charged particle 
oscillates near the upper conductor plate in Fig. [2] 



F. Electromagnetic radiation 

It is well known that the oscillating charged-particle radi- 
ates electromagnetic energy. With respect to the reference 
point on the surface of the upper conductor plate, the oscil- 
lating charged-particle has a dipole moment given by 
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and 



Pd = -eAneoK^vs 



where = (t) , s = s (t) , and Eq. d84l > has been inserted for 
Qj. The negative sign comes from the fact that the particle is 
perceived as residing in the negative z-axis to someone on the 
surface of the upper conductor plate. In terms of z c i = Zd (?) 
defined in Eq. d90l ), p c / becomes 

Pd = -e z 47teoK 3 v(zd+b), 

Pd = -e z 47Z£ K 3 VZd, 

where b is a constant. The electromagnetic power radiated by 
an oscillating charged-particle, P rac i, is given by the Lienard 
formula, 
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Insertion of Eq. d9TT l for 'id finally yields 

87T£oK- 2 V 2 / Treo^v J v 



Prad 



3c 3 ^ 4m \(z d + b) 2 
V [y(b 3 -a 3 )-b 3 ]E p 



(h-zd-by 



(zd + bY 



where 



(h-Zd-by 
2a{b~a)G\ a 2 0\+b 2 02 



eoK 2 £qK 3 

which is the result defined in Eq. (fT7b for the Lienard radiation 
power. 

IV. CONCLUDING REMARKS 

The phenomenon of charged-particle oscillation subjected 
to a constant electric field has been investigated. For a posi- 
tively charged core-shell structured particle, the criterion for 
an oscillatory motion is given by 
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For the parallel plate configuration illustrated in Fig. [2] the 
positively charged core-shell structured particle can only have 
oscillatory modes in region where < Zd < h/2. This region is 
further divided into subregions A and B, and this is illustrated 
in Fig. [8] In this work, the upper conductor plate is located 
at Zd = and the lower conductor plate is located at Zd = h. 
Further, the upper conductor plate is held at voltage Vj and 
the lower conductor plate has a voltage of Vl, where Vj >Vl- 
That said, the magnitude of a dominant force in region A falls 
off with distance like 1/s 3 and this force acts to repulse the 
positively charged core-shell structured particle from the up- 
per conductor plate. In region B, the magnitude of a dominant 
force falls of like 1/s 2 with distance and this force acts to at- 
tract the particle towards the upper conductor plate. It is this 
competition between the two dominant forces from regions A 
and B that puts particle in an oscillatory motion. 

Such oscillatory mechanism does not involve any charge 
exchange; therefore, it is fundamentally different from the tra- 
ditional description by charge exchange. 9 Nevertheless, the 
novel finding in this work does not invalidate the charge ex- 
change description altogether because the two oscillatory phe- 
nomena are not exactly the same. For instance, the traditional 
picture deals with charged particle oscillation in which the 
particle sweeps the entire gap between the two plates whereas, 
here, the particle only oscillates in the restricted regions be- 
tween the plates. Further, here, the charged-particle must be 
structured and not a "point" particle in order to have any os- 
cillatory motion. The reason for this is because the repulsive 
force in region A arises as a consequence of induced polariza- 
tion and a point particle does not have such property. Without 
such force appearing in region A, there would not be a way 
to repulse the positively charged particle from sticking to the 
surface of the upper conductor plate in Fig. [8] The traditional 
description by process of charge transfer does not have such 
restrictions, of course. Thus, the two phenomena are not ex- 
actly the same. Nonetheless, the novel finding here presents 
yet another mechanism for a charged-particle oscillation in an 
otherwise uniform and constant electric field. 

For a negatively charged core-shell structured particle, the 
oscillatory criterion is given by 
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where 
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The oscillatory behavior for negatively charged core-shell 
structured particle is observed only in the region where h/2 < 
Zd < h; and, this is illustrated in Fig. [12] For the case of neg- 
atively charged particle, subregions A and B are formed near 
the lower conductor plate side of Fig. Q~2j which is exactly 
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the opposite of the positively charged particle case. Again, 
the force in region A repulses the particle from the lower con- 
ductor plate and the force in region B attracts the particle to 
the lower conductor plate. It is this "push-pull" competition 
between the two dominant forces in regions A and B that puts 
charged-particle in an oscillatory motion. Such oscillatory be- 
havior is only possible because the particle can be polarized 
under applied electric field. 

Because of the explicit mass dependence in Eqs. ( llOOl i and 
( llOU , the charged particle oscillation discussed in this work 
is more favorably satisfied by microscopic or smaller parti- 
cles than by macroscopic counterparts due to various experi- 
mental limitations. In principle, particles of any size can be 
charged to satisfy the oscillation criterion of Eqs. (1 100b or 
(1101b provided a DC electric field of sufficient strength can 
be applied without electrical breakdown taking place. Unfor- 
tunately, electrical breakdown occurs at some point even in 
vacuum and this may limit macroscopic particles from satis- 
fying Eqs. ( fTOOl or (ITOB^- 

The charged-particle oscillator based on the presented 
novel mechanism represents a natural prototype for illuminat- 
ing electric dipole radiation. In such system, the frequency 
of emitted electromagnetic radiation is controlled by a DC 
voltage biased across the two plane-parallel electrodes. The 
strength of emitted radiation power from such system depends 
directly on the magnitude of effective charge carried by the 
charged-particle. 

As for potential applications, the finding in this work can 
be utilized to build a source for generating microwave radia- 
tion. Microwaves thus generated, for instance, might be used 
to heat water or to excite gases in tiny plasma capsules to 
produce light. As another potential application, the repulsive 
mechanism discussed here can be utilized to build an anti- 
friction device. Such device would not require any grease or 
oil, which are known to be environmentally hazardous. In- 
stead, the friction in such device would be controlled by ap- 
plied DC voltage. 

In nanoscale plasma confinements, nanocavities for in- 
stance, only scarce number of ionized atoms or particles par- 



ticipate in the dynamics. This is so because only finite number 
of ionized particles can be contained in such small scale con- 
finement. In such system, effects arising from the collective 
motion of particles may only be a small perturbation com- 
pared to the effects arising from individual particle dynam- 
ics. For very large scale systems, wherein enormous num- 
ber of ionized particles are involved, exactly the opposite is 
true. There, effects arising from the collective motion of 
particles become predominant. As the technology evolves, 
plasma systems will eventually enter the regime of submicron 
to nanometer scale plasma confinements. Already, the size 
of display pixels is approaching the dimensions of few visi- 
ble wavelengths; and, display technologies based on charged 
particles are beginning to involve fewer ionized nanoparticles, 
where individual particle effects are no longer small perturba- 
tions. In the electronics industry, the finite particle nonneutral 
plasma systems are becoming technologically very important 
as nanoscale fabrication processes demand for the develop- 
ment of extreme UV (ultraviolet) to X-ray lasers. For in- 
stance, device fabrication process at length scale of ~ 20 nm 
or less by photolithography requires extreme UV lasers. Be- 
cause solid state devices cannot generate such high frequency 
laser waves, plasma based sources are the only viable candi- 
dates for building extreme UV or X-ray lasers. The macro- 
scopic plasma sources can readily generate extreme UV elec- 
tromagnetic waves, but at the cost of losing the coherent na- 
ture of laser waves. Inevitably, the coherent nature of a laser 
source demands for finite particle plasma sources, wherein 
the physics of individual particle dynamics is predominantly 
important. In this respect, the results obtained in this work 
should be useful and interesting to certain areas of plasma 
physics. 
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